ON THE RIGOROUS DERIVATION OF THE 3D CUBIC NONLINEAR 
SCHRODINGER EQUATION WITH A QUADRATIC TRAP 
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Abstract. We consider the dynamics of the 3D A'^-body Schrodinger equation in the presence of a quadratic 
trap. We assume the pair interaction potential is N^^~^V{N^x). We justify the mean-field approximation 
and offer a rigorous derivation of the 3D cubic NLS with a quadratic trap. We establish the space-time 
bound conjectured by Klainerman and Machedon [30] for /3 G (0,2/7] by adapting and simplifying an 
argument in Chen and Pavlovic [7] which solves the problem for /3 G (0, 1/4) in the absence of a trap. 
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1. Introduction 

It is widely believed that the cubic nonlinear Schrodinger equation (NLS) 
O' idt(f> = L(l) + \(l>\^ (f), 



(N 



X 



where L is the Laplacian — A or the Hermite operator — A + w'^ |a;|^ , describes the physical phenomenon of 
Bose-Einstein condensation (BEC). This belief is one of the main motivations for studying the cubic NLS. 
BEC is the phenomenon that particles of integer spin ("Bosons") occupy a macroscopic quantum state. This 
unusual state of matter was first predicted theoretically by Einstein for non-interacting particles. The first 



^ I experimental observation of BEC in an interacting atomic gas did not occur until 1995 using laser cooling 
■ " " techniques [1, 14]. E. A. Cornell, W. Ketterle, and C. E. Wieman were awarded the 2001 Nobel Prize in 

Physics for observing BEC. IMany similar successful experiments were performed later on [4, 13, 27, 37]. 

Let t G R be the time variable and xjv = {xi, X2, xn) G ^'^^ be the position vector of N particles in 
M.^. Then BEC naively means that the N-body wave function ipj^{t,'XN) satisfies 

N 

V'A,(i,XA,) = J]^ 0(i,a;j) (1.1) 

up to a phase factor solely depending on t, for some one particle state (p. In other words, every particle is in 
the same quantum state. It is believed that the one particle wave function (/) which models the condensate 
satisfies the cubic NLS. Gross [25, 26] and Pitaevskii [34] proposed such a description. However, the cubic 
NLS is a phenomenological mean field type equation and its validity needs to be established rigorously from 
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the many body system which it is supposed to characterize. As a result, we investigate the procedure of 
laboratory experiments of BEC according to [1, 14]. 

Step A. Confine a large number of Bosons inside a trap e.g., the magnetic fields in [1, 14]. Cool it down so 
that the many body system reaches its ground state. It is expected that this ground state is a BEC 
state / factorized state. This step corresponds to the mathematical problem. 

Problem 1. Show that the ground state of the N-hody Hamiltonian 

j=l ^ ' l^i<j^iV 

is a factorized state. 

2 

We use the quadratic potential \x\ to represent the trap. This simplified yet reasonably general 
model is expected to capture the salient features of the actual trap: on the one hand the quadratic 
potential varies slowly, on the other hand it tends to oo as — > oo. In the physics literature, 
Lieb, Seiringer and Yngvason remarked in [31] that the confining potential is typically ~ |x| in the 
available experiments. Mathematically speaking, the strongest trap we can deal with in the usual 
regularity setting of NLS is the quadratic trap since the work [39] by Yajima and Zhang points out 
that the ordinary Strichartz estimates start to fail as the trap exceeds quadratic. 
Step B. Switch the trap in order to enable measurement or direct observation. It is assumed that such a 
shift of the confining potential is instant and does not destroy the BEC obtained from Step A. To 
be more precise about the word "switch": in [1, 14] the trap is removed, in [37] the initial magnetic 
trap is switched to an optical trap, in [4] the trap is enhanced, in [13] the trap is turned off in 2 
spatial directions to generate a 2D Bose gas. Hence we have a different trap after the switch, in 
other words, the trapping potential becomes The system is then time dependent unless 

LJ = ujQ- Therefore, the factorized structure obtained in Step A must be preserved in time for the 
observation of BEC. Mathematically, this step stands for the following problem. 

Problem 2. Take the BEC state obtained in Step A. as initial datum, show that the solution to 
the many body Schrddinger equation 

ia,V;v = E(-^A.,+^'^)V';v + ^ E N'^V{N^{xi-xj))i^^ (1.2) 

j = l \ / l^i<j!iN 

is a BEC state / factorized state. 

We first remark that neither of the problems listed above admits a factorized state solution. It is also 
unrealistic to solve the equations in Problems 1 and 2 for large N. Moreover, both problems are linear 

so that it is not clear how the cubic NLS arises from cither problem. Therefore, in order to justify the 
statement that the cubic NLS depicts BEC, we have to show mathematically that, in an appropriate sense, 

N 

'iljj^{t,'x.N) JJ (j){t,Xj) as N ^ oo 
for some one particle state (j) which solves a cubic NLS. However, when 4> 4>' 

2 

N N 

Y\_<p(t,Xj) - Y\_'p(fi^j) -^2asN^oo. 
i.e. our desired limit (the BEC state) is not stable against small perturbations. One way to circumvent this 

(k) 

difficulty is to use the concept of the k-particle marginal density j)^' associated with ip^ defined as 

7iv^(i,Xfc;xJ.) = j Viv(*)Xfc,xjv-fe)'^i'jv(i.Xfe,XAr_fe)dxjv-fe, x^jX^eM^'' (1.3) 

and show that 

k 

7j^^(t,Xfc;x'j,) ~ Yl(t>{'t,Xj)^{t,x'j) as A/' oo. 

i=i 

Penrose and Onsager [33] suggested such a formulation. Another approach is to add a second order correction 
to the mean field approximation. See [11, 23, 24]. 

For Problem 1, Lieb, Seiringer, Solovej and Yngvason showed that the ground state of the Hamiltonian 
exhibits complete BEC in [32], provided that the trapping potential Vtrapix) satisfies inf Vtrap{x) — >■ oo 
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for i? — oo and the interaction potential is spherically symmetric. To be more precise, let V'jvo t>e the 
ground state, then 

^n]o \<t>Gp) {<t>Gp\ as AT ^ DO, 

where is the corresponding one particle marginal density defined via formula 1.3 and (f)Qp is the 

minimizer of the Gross-Pitaevskii energy functional 



/ 



( I V(/.|' + Vtrap{x) |</.|' + 47rao l-^l^ )dx. 



So far, there has not been any work regarding the Hamiltonian evolution in Step B in the case when 
w ^ 0. Motivated by the above considerations, we aim to investigate the evolution of a many-body Boson 
system with a quadratic trap and study the dynamics after the switch of the trap. We derive rigorously the 
3D cubic NLS with a quadratic trap from the AT— body linear equation 1.2. To be specific, we establish the 
following theorem in this paper. 

Theorem 1. (Main Theorem) Let \i^n^ family of marginal densities associated with ipj^, the 

solution of the N-body Schrddinger equation 1.2 for some (3 G (O, |] . Assume that the pair interaction V is 
a nonnegative L^{M.^) n n W^'°°{M.^) spherically symmetric function. Moreover, suppose the initial 

datum of equation 1.2 verifies the following the conditions: 
(a) the initial datum is normalized i.e. 

IIV'ivWIL. = 1, 

(h) the initial datum is asymptotically factorized 

lim Tr 7^^(0,a;i;a;i) - 0o(a;i)0o(a;i) =0, 

A/— 7- CO 

for some one particle wave function (pQ . 

(c) the initial datum has hounded energy per particle i.e. 

sup ^ (V'jv(O), HniPn{^)) < oo> 

where the Hamiltonian is 

j=l \ / l^Kj^N 

Then Vi ^ 0, Vfc > 1, we have the convergence in the trace norm that 

k 



lim Tr 

iV->cx> 



7j^^(i,Xfe;x'fc) - Yl 4>{t,Xj)(j){t,x'^ 

where (p(t, x) is the solution to the 3D cubic NLS with a quadratic trap 

idt<i> = ||-^A, + w2^j </) + 6o|<^|''/' mK3+i (^ ^^ 

(f){Q,x) = <i)Q{x) 
and the coupling constant bo = J^a V{x)dx. 

For the u) = case, the approach which uses the marginal densities |77v'} dynamics problem 

has been proven to be successful in the fundamental papers [15, 16, 17, 18, 19, 20, 21] by Elgart, Erdos, 
Schlein, and Yau. As pointed out in [21], their work corresponds to the evolution after the removal of the 
traps. Motivated by a kinetic formulation of Spohn [36], their program consists of two principal parts: in 

one part, they prove that an appropriate limit of the sequence ji^' | as A'' ^ 00 solves the Gross-Pitaevskii 
hierarchy 

(iat + iA.,-iA.,)7^^^=«'oE^^-.fe+i(^^'^'0' ^ = 1'-'"'- (1-5) 
where B^^k+i = B]^^^^ - B'l^^^^ and 

B],k+i (7^''"^^^) (i,Xfc;x^) = j j - Xk+i)5{xk+i - x'k+i)-i^^^^\t,y.k+i;'x'k+i)dxk+idx'k+i, 
Bl^+^ (7^''+'^) (i,Xfc;x;.) = j y"(5(x;-Xfc+i)5(xfc+i -4+i)7('=+i)(t,Xfe+i;x'fe+i)dxfe+id4+i; 



4 



XUWEN CHEN 



in another part, they show that hierarchy 1.5 has a unique solution which is therefore a completely factorized 
state. However, as remarked by Terence Tao, the uniqueness theory for hierarchy 1.5 is surprisingly delicate 
due to the fact that it is a system of infinitely many coupled equations over an unbounded number of 
variables. In [30], by imposing a space-time bound on the limit of Ittv^I' Klainerman and Machcdon 
gave another proof of the uniqueness in [18] through a collapsing estimate originated from the ordinary 
multilinear Strichartz estimates in their null form paper [29] and a board game argument inspired by the 
Feynman graph argument in [18]. 

Later, the method in Klainerman and Machedon [30] was taken up by Kirkpatrick, Schlein, and StafHlani 
[28], who studied the corresponding problem in 2D; by Chen and Pavlovic [5, 6], who considered the ID 
and 2D 3-body interaction problem and the general existence theory of hierarchy 1.5; and by the author 
[12], who investigated the trapping problem in 2D. In [8, 9], Chen, Pavlovic and Tzirakis worked out the 
virial and Morawetz identities for hierarchy 1.5. 

Recently, for the 3D case without traps, Chen and Pavlovic [7] proved that, for (3 € (0, 1/4), the limit 

of |77v^| actually satisfies the space-time bound assumed by Klainerman and Machcdon [30] as iV — > 
00. This has been a well-known open problem in the field. Moreover, they showed that the solution 
to the BBGKY hierarchy converges strongly to the solution to hierarchy 1.5 in without assuming 
asymptotically factorized initial datum. In this paper, we adapt and simplify their argument in establishing 
the Klainerman-Machedon space-time bound. We also extend the range of (3 from (0, 1/4) in Chen and 
Pavlovic [7] to (0, 2/7] . Through simple functional analysis, we obtain a convergence result without assuming 
asymptotically factorized initial datum as well. (See Corollary 1) But we are not claiming it as a main result 
in this paper. We compare our result and the one in Chen and Pavlovic [7] briefiy in Section 1.1 below. 

1.1. Comparison with Chen and Pavlovic [7]. For comparison purpose, we transform Theorem 4 in 
Section 2, which implies Theorem 1, into the general convergence result below without the assumption of 
asymptotically factorized initial data since the result in Chen and Pavlovic [7] is under the same regularity 
setting (condition 1.6) as Theorem 4. 

Corollary 1. (For comparison purpose only, not a main result.) Let 7^'' (t, x^; xj.) he the k-marginal density 
associated with the solution of the N-body Schrddinger equation 1.2 for some (3 G (0,2/7] and "i'^''\t,Xk'.x'i^) 
be the solution to the Gross-Pitaevskii hierarchy with a quadratic trap (hierarchy 3.3). Assume the initial 
data of and 7^*^^ satisfy the following conditions 
(a') 



lim Tr 



7iv'^(0,Xfc;x',)-7('='(0,Xfc;x',) 



(b'J 



Then we have the convergence of the evolution in trace norm 



(1.6) 



lim Tr 

N->-oo 



7iv'\i,Xfe;x',)-7('=)(t,Xfc;x',; 



= 0. 



Proof. This result for the non-trap case should be credited to Erdos, Schlein and Yau though they did not 
state it as one of their main theorems. They mentioned it on page 297 of their paper [21]. Once we have 
established Theorem 4, its proof together with some simple functional analysis proves this corollary. We 
include the proof in Appendix I (Section 9) for completeness. □ 



Briefly, the main result (Theorem 3.1) in Chen and Pavlovic [7] is the following. 

Theorem 2. [7] Let u = 0. Suppose 7^'' (f, x^; x^) is the k-marginal density associated with the solution of 
the N-body Schrddinger equation 1.2 for some (3 G (0, 1/4) and 7*^*^-' (t, x^; x'^.) is the solution to the ordinary 

Gross-Pitaevskii hierarchy without traps. If the initial data of 7^^ and 7^*^^ satisfy the following conditions 
(a") 



lim 

N—KX 



n (1 - A.J*+^ (1 - A.,)"^' (7iv'^(0,x,;x',) - 7('=)(0,x,;x',) 



= 0, 



L^(dx.kdx'.) 



(b") 



<V;v(o),<V'iv(o))^c'=iv^ 
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then we have the convergence of the evolution in norm 



lim 

N-^oo 



= 0. 



One can then easily tell the following: on the one hand, the result in this work allows a quadratic trap 
and a larger range of /3 in the analysis (these are also the main novelty and the main technical improvement 
of this paper); on the other hand, the result in [7] yields a stronger convergence {H^ convergence) when 
the initial data admits a stronger convergence. The main purpose of this paper is to justify the mean-field 
approximation and offer a rigorous derivation of the 3D cubic NLS with a quadratic trap (Theorem 1). Thus 
we establish the convergence of probability densities which is the trace norm convergence in this context. 
The Chcn-Pavlovic result is crucial for their work on the Cauchy problem of Gross-Pitaevskii hierarchies. 
Before we delve into the proofs, we remark that taking the coupling level k to be IniV in Section 6.0.3 is 
exactly the place where we follow the original idea of Chen and Pavlovic [7]. 

1.2. The Anisotropic Version of the Main Theorem. It is of interest to use anisotropic traps in 
laboratory experiments. (Sec, for example, [13].) Our proof for Theorem 1 also applies to the case with 
anisotropic traps. In fact, wc have the following theorem. 

Theorem 3. Define the diagonal matrix 



Q 



Let Itjv"*} the family of marginal densities associated with ip^^, the solution of the N-body Schrodinger 



equation 

^ / 1 1 \ 1 

idti^N = E - 2^-0 + Q^j + TV 2Z N^^V{N^ (ar, - ar,))V^ 
i=i ^ ^ iii%<j^N 

for some /? G (O, |]. Assume that the pair interaction V is a nonnegative L^{R^) n H^{M.^) n W^-'°^( 
spherically symmetric function. Moreover, suppose the initial datum of equation 1.2 verifies Conditions (a), 
(b), and (c) in the assumption of Theorem 1. Then Vt ^ 0, VA; ^ 1, we have the convergence in the trace 
norm that 

k 

hm Tr ^f{t, Xfe; x^) - [] (^{t, Xj)4>{t, x'^) = 0, 
where ^{t, x) is the solution to the 3D cubic NLS with anisotropic quadratic traps 



idt4> 



1 

' L 

2 



1 



Qx](j) + bo\<f)f<p mR^+^ 



and bo = /^a V{x)dx. 

The anisotropic version of the main theorem stated above yields to the same techniques as Theorem 1 
but is technically more complicated. Therefore we prove only the latter in detail, merely suggesting during 
the course of the proof the appropriate modifications needed to obtain the more general theorem. 

2. Proof of the Main Theorem (Theorem 1) 
We establish Theorem 1 with a smooth approximation argument and the following theorem. 

Theorem 4. Let |7jv'} family of marginal densities associated with ipj^, the solution of the N- 

body Schrodinger equation 1.2 for some e (0, f ] . As sume that the pair interaction V is a nonnegative 
I/^(M'^) n if^(]R^) n W'^'°°{E?) spherically symmetric function. Moreover, suppose the initial datum of 
equation 1.2 is normalized, asymptotically factorized and verifies the condition that there is aC independent 
of k or N such that 

(Vjv(o),irW;v(o))^c'=iv'= 

Then Vt > and Vfc > 1, we have the convergence in the trace norm that 



lim Tr 

N—^oo 



7^^(i,Xfe;x'fc) - Yl4){t,Xj)(j){t,x'j] 
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where 4>{t, x) is the solution to the 3D cubic NLS with a quadratic trap 

I |2 \ 



idt4> 



m 



and bo = J^s V{x)dx. 

Before presenting the proof of Theorem 4, we discuss how to deduce Theorem 1 from Theorem 4. It is 
a well-known smooth approximation argument. We include it for completeness. For technical details, we 
refer the readers to Erdos-Schlein-Yau [20, 21] and Kirkpatrick-Schlein-Staffilani [28]. 

Write the spaces of compact operators and trace class operators of L^(R'^'^) as K-k and Then (/Cfe)' = 
Cl- Via the fact that /C^ is separable, we select a dense countable subset of the unit ball of /Cfe and call it 



C /Cfc. We have 



4, for7W,7('=)e4: define 



< 1 where ||-|| is the operator norm. We set up the following metric on 



E2- 

i=l 



TtJ 



Then a uniformly bounded sequence € £^ converges to 7^'^^ € with respect to the weak* topology 
if and only if 

lim dfe(7iv'\7''^) = 0. 

N^oc 

Fix K > and x e W, with < x < 1, x(s) = 1, for ^ s < 1, and x(s) = 0, for s > 2. We 
regularize the initial data of the A/'-body Schrodinger equation 1.2 with 

^^^^ ||x(«iJiv/A^) Vat (0)11 ' 
and we denote ipjs[{t) the solution of the A^-body Schrodinger equation 1.2 subject to this regularized initial 
data and < j)^ (t) > the family of marginal densities associated with 'ipj^{t). 

L J k—l 

With these notations, if k > small enough, on the one hand, we have 

(V-AT (0),if^V'Ar (0)) ^(7'=iV^ 

and 



lim^TrjW |7^^)(O)-n'^o0ol = 0, 



for every J^'^^ £ K-k- Thus 7!^' (0, x^; xj.) TT 4>a4'o as ^ 00 in the weak* topology. Since TT </'o'?^o 

is an orthogonal projection, the convergence in the weak* topology is equivalent to the convergence in the 
trace norm. Consequently, the conditions of Theorem 4 are verified and it implies that Vt e [0, Tq] and 
Vfc > 1, 

k 

7i^^(*,Xfc;x'fc) - Y[ (l>it,Xj)'^{t,Xj] 

i=i 

On the other hand, there is a constant C independent of N and k such that 



lim Tr 



= 0. 



Tr 



J''''{l'^\t)-j^^\t)) 



for every Jf'^) e /Cfe. Therefore, 



Tr J« j'-^\t,^k;^',) -l[cPit,Xj)c^{t,x'j) 



Tr 



^ C 



j(fc) 



+ 



TrjW [j'-;^\t)-l[ct>{t,xj)cl>{t,x'^) 



+ 



TrjC'^ [j'-^\t)-l[(t>{t,xj)<t>{t,x'j) 
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As the above inequality holds for all k > small enough, wc know 7l^-'(t. xj,; x',,) TT (l){t,Xj)(p{t.x'A 

as A'^ — > oo in the weak* topology. This convergence is again equivalent to the convergence in the trace 

norm because I I <f){t,Xj)(j){t,x'A is an orthogonal projection as well. Whence we have established our 

main theorem (Theorem 1) through Theorem 4. It remains to show Theorem 4. We prove Theorem 4 with 
the help of the lens transform. We define the lens transform and show its related properties in the next 
section then we establish Theorem 4 in Section 4. 

3. Lens Transform 

In this section, we first define the lens transform and review its relevant properties, then we prove an 
energy estimate (Proposition 2) which relates the energy on the two sides of the lens transform. It aids 
in the proof of Theorem 4 in the sense that it links the analysis of — A^, + uP' to the analysis of —Ay 
which is a better understood operator. We denote (f, x) the space-time on the Hermite side and (r, y) the 
space-time on the Laplacian side. We now define the lens transform we need. 

Definition 1. Let 'x.n,Yn € W^^ . We define the lens transform for L"^ functions Mm ■ L^{dyN) — > 
L^{dxN) and its inverse by 

I 1^ 

(MjvUjv) (t,Xjv) = — — 7j — un{ , -) 

(cosojt)— w coscjt' 

1 , \ / ^ e n-"^!-^ 2 arctan(a;r) yn 
M^Viv (r,yiv) = —-^i'Ni -^,^=^=). 

Mjv is unitary by definition and the variables are related by 

tanwf Xfc 



LO COS Ldt 

Definition 2. Let Xfe,x^,yfe,y^ e M"*^. We define the lens transform for Hilbert- Schmidt kernels : 
L^idykdy'i-) L'^{dxkdx'i.) and its inverse by 

v / (coswt)"* u) coscot coswr 

. ^2 (|yfc|^-KP) 

^ (1-1- oj^r^^j 2 w vl + 1^ T-' vl + w^T^ 

Tfc is unitary by definition as well and the variables are again related by 

tanwi Xfc , , x'. 
T = , Yfc = 7 and = 



U! COS Ojt COS Ljt 

Before we characterize the accurate effect of the lens transform, we clarify the motivation of such defini- 
tions by a lemma. 

Lemma 1. [3, 12] Define a and (5 via the system 

a{t) + r]{t)a{t) = 0, a(0) = 0, a{0) = 1, 

m + r,{t)m = 0,/3(0) = l,/3(0)=0. 

If /3 is nonzero in the time interval [0,T], then the solution of the ID Schrddinger equation with a time 
dependent quadratic trap 

idti) = (^-^dl + ^T]{t)x^^ in 
ilj{0,x) = f{x)€L\R) 

in [0, T] is given by 

e flW 2 r<e(t) X , 
tp{t,x) = rW(-57TT, ^TTT , 



if u{t, y) solves the ID free Schrddinger equation 
subject to the same initial data. 



id^u = -^d^u in ]R^+^ 
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Proof. See [3] for a proof by direct computation and [12] for an algebraic proof using the metaplectic 
representation. When r]{t) = w^, such a transformation has a long history, we refer the readers to [3] and 
the references within. □ 

To make formulas shorter, let us write 

(t,-SW) (r,yfe;3^) =7^'=)(i,Xfc;x',)/iW (r,y,;y;,) 

then precisely, the lens transform has the following effect. 

Proposition 1. 

hl^Hr,yk;y'k) f.r, ( 1, ,_o|xfe|^\ / 1 ,,,2l4l^ 



(1+WV2) 

Proof. This is a direct computation. □ 

Via Proposition 1, wc know how the lens transform acts on the Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy and the Gross-Pitaevskii hierarchy. 

Lemma 2. {BBGKY hierarchy under the lens transform) Write Vn {x) = N^^V (N^x) . |7jv^} ^'^^'ves the 
BBGKY hierarchy with a quadratic trap 

■a (k) I Ia , 2|Xfc|^\ (fc) , / 1a , 2|Xfcl^ 1 (fe) /oi\ 

^daV - I - 2 Axfc + w ^- jTN +\ - + 2 J ^ 

H — ]y~E / (^Jv(a;j - Xfe+i) - Viv(a;j- - Xfe+i)) 7^^"^^^ {t,y:k,Xk+i;x'k,Xk+i) dxk+i, 

m [—To, To] if and only if ^u^^^ = Tj^^7j^''| solves the hierarchy 

iSr + ^Ay.-iAy,)^) (3.2) 



in 



r tan coTq tan ojTq 1 



Lemma 3. {Gross-Pitaevskii hierarchy under the lens transform) {7^'^^} solves the Gross-Pitaevskii hier- 
archy with a quadratic trap 



in [—To, To] if and only if {u*^'^^ = Tj^^7('')} solves the hierarchy 

(idr + ^Ay, I Ay,] u^'^) = L__6o^B._,+in(*=+i), (3.4) 



in 



r tan wTn tan cuTp 1 

L w ' w J * 



The lens transform for the Hilbcrt-Schmidt kernels is not only by definition a unitary transform on 
(]^2nfe^ ^ g^^gQ g^j^ isometry on the space of the trace class operator kernels. 
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Lemma 4. £ [— ro,To], Viir(y/c,y^) the kernel of a trace class operator on (M"*^). // 
then 

' {TkK) (xfe, x'fe) {MmJ) (xi) dx'fe = A (M^/) (x^) . 



/n of/ier words, the eigenvectors of the kernel (Tfei^) (xfe,xj^) are exactly the lens transform (Lemma 1) of 
the eigenvectors of the kernel K{yk,y'i.) with the same eigenvalues. In particular, we have 

Tr \TkK\ =Tr\K\. 

Proof. This is a straight forward computation. Wc remark that we have defined the generaUzed lens 
transform for a function and a kernel separately via Definitions 1 and 2. □ 

Once we have proved the following proposition which relates the energy of the two sides of the lens 
transform, we can start the proof of Theorem 4. 

Proposition 2. Let '4'jq{t, xat) he the solution to equation 1.2 for some j5 € (0, 3/5) subject to initial '4'n{^) 
which satisfies the energy condition 

If UN{T,yN) = Afjv^V'AT; then there is a C ^ 0, for all fc > 0, 3A^o {k) such that 

{ UN 

for all N^No and all r e [-isji^, tan^j provided that Tq < ^. 

The rest of this section is the proof of Proposition 2. We prove it for w > through Lemmas 5 and 6 
since the case to = has already been studied in [15, 18]. 

Lemma 5. For (3 G (0, 3/5) , there is a C ^ 0, for all k ^ 0, 3No (fc) such that 



{^,H%^) > C'^N^ /<p,I] (-A,. ) 



for all N ^ No and all (p & L^ 



Proof. The proof basically follows Proposition 1 in [15] step by step if one replaces (l — A^;^) by ^— 
and notices that for w > 0, 



+ oj^ \x 



(1-A)^/ 



(- 



([38]) 



when one uses Sobolev. There are some extra error terms which can be easily handled. We illustrate the 
control of the extra error terms through the following example. Write Sf = ^— A^,. + consider 

N^^ {{^, S!...Sl+,V {N^ (xi - X2)) + cc) 
= TV^^ {{S^...Sn+w, SlSlV [N^ {xi - X2)) S3...Sn+w) + cc.) , 
where cc. denotes complex conjugates. Neglecting Ss...Sn+i, we have 

N^P {{^, SfSiV {N^ (xi - X2)) ^) + cc) 
= N^^ ((^, (-A,, + uj^ laril") (-A,, + oj^ \x2f) V {N^ (an - ^2)) + cc.) 
= N^^ {(ip, (-A,J (-A,J V [N^ {xi - X2)) if) + cc) 

+N^^ ((^, (-A,J (oj^ \x2f) V {N^ {xi - X2)) if) + cc) 

+N^^ ((<p, (oj^ \xi\') (-A,,) V (N^ (xi - X2)) if) + cc) 

= I + II + III + IV 
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Compared to [15], the extra error terms are //, ///, and IV. Since IV is positive, we only look at 

II ^ -2N^^\{V^,Lj\x2\ip,\{VV){N^ {xi-X2))\L0\x2\ip)\ 

> -CN^I'a (V,,w \X2\ I (Vy) {N^ {xi - x^)) \ V ^.u \x2\ if) 
-CN^^a-^ {to \X2\ ip, I (Vy) {N^ (xi - X2)) \uj\x2\ 

> -CiV4/3 (a (v3, SlSlip) + a-'^N-'^l^ (w \x2\ S'^w |x2| ^)) (Sobolev at the second term) 

> -CN''^{^, SlSlip) {a = N-^) 

^ -CN^^-^N^{^,SfSl^) 

As long as ^ e (O, |) , we can absorb the extra error terms into the main term N'^ (^ip, SfS^f)- The Sobolev 
we used is 



J V {N^ {xi - X2)) dxidx2 < \\V{N'^-) 



3 



\ip\ dxi I dx2 



L2 



dxidx2 



2 I |2\ 2 



dx\dx2- 



□ 



Lemma 6. iet 



-fx (i) = i^x cos cot — cox sin Lot. 



//ujv(r,yjv) = (i/;jv), t/ie^ 



{un{t), {-Ay.) un{t)) = (P,, (i) (t) V-^W) 



Proof. We provide a proof through direct computation here. We remark that (t) is in fact the evolution 
of momentum. See [12]. Without lose of generality, we may assume N ~ 1, then 



P^(i)Vi(i) = Px(i) 



e "2°" '^1 / tanwt 



{cos Loty \ ^ [cos ujt) 



(coswt) ^ 

e 2 



(iVx COSWt) Ml 



tan a; 



LJ ' {cos LOt) 



(cos wt) ■ 



(zV3^)wi (t,2/), 



Thus 



(t.i(T),(-A,)ni(r))= / \VyU^{T,y)fdy= [ ^^'^^'' e^or = (P. ft) Vift),Px ft) ViW) • 

J 7 coswt 



□ 



With the above lemmas, we prove Proposition 2. We first notice that 



un{t), n (1 - A^.) UN{r)\ = / Vivft), n (1 + ft)) ^ivft) ) (Lemma 6) . 



Since 



= (/) /) + cos^ Lot (V/, V/) + sin^ {xf, xf) — 2io sin cot cos Im (V/, xf) 
< c(/,(-A,+a;2lxl2)/) 



we know 



V'ivft),n(i+<- ^ c'=(v'ivft),n(-^-.+^'i^^i')^^w 



^ {^N{t),H^i^Nit)) ■ (Lemma 5) 
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From the energy condition on the initial datum tp]^{0), we then deduce Proposition 2 that is 

( Un{T), 

4. Proof of Theorem 4 

We devote this section to estabUshing Theorem 4. The main idea is to first prove that, in the time 

period r e [O, i^iL^] ^ith To < ^, as N oo, = T^"Si^^} , the lens transform of the solution to 

the BBGKY hierarchy 3.1, converges to {u^'^' = Tj^^'y^''^}, the lens transform of the solution to the Gross- 

Pitaevskii hierarchy 3.3 in the trace norm, then use Lemma 4 to conclude the convergence 7^^' 7''^^ in 
the trace norm as -/V 00 in the time period t G [0,To]. Then the time translation invariance of equation 
1.2 proves Theorem 4 for all time. We now present the proof in detail. 

Step I. In Proposition 2, we have already established the energy estimate for {ujv = M^^tp^} , 

sup (uNiT),Y[{i-^y,)uN{T)) 
re[0,'-^] \ j = l / 

which becomes 

sup Trlr\{l-Ay^)]u^^^^C\ (4.1) 

for |uj^^ = T/^^^N^^- Therefore we can utilize the proof in Erdos-Schlein-Yau [18] or Kirkpatrick- 
Schlein-Staffilani [28] to show that the sequence jw^' = T^^^^^ | is compact with respect to the 
weak* topology on the trace class operators and every limit point {u^*^^} solves hierarchy 3.4 which 



IS 



(idr + ^Ay, - ^Ay, ) = yrT^6o Bj,k+iu^>^+'\ (4.2) 

in 3D. This is a fixed time argument. We omit the details here. 
Notation 1. To make formulas shorter, let us write 

5(r) = (l + c.V)-^ 

from here on. The only property of g{r) we are going to need is that < c < g{T) < C < 00 m any finite 
time period. 

Step II. In this step, we use the a-priori estimate 4.1 to provide a space-time bound of u^''^ so that we can 
employ Theorem 6 in Step III. We transform estimate 4.1 into the following theorem. 

Theorem 5. (Main Auxiliary Theorem) Assume V is a nonnegative L^(M^) n iJ^(M^) function. Let 

VnAv) = N^^Vr{N^y) = N^^g^{T)V{g{T)N^y) 

be the interaction potential with the interaction parameter /3 e (O, |] . Suppose that uj^^ solves hierarchy 3.2 
in [0, T] c [0, ifiJL^iZii] ^ which, written in the integral form, is 

•^'^ = (r)4?o (4.3) 



1 -^0 y\ J 



k 



N-k 



N 

3 



subject to the condition that 



sup Tr \\{{l-Ay^) uP^C\ (4.4) 



re[0,T] 
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where BN,j,k+i,T = -BArj,fc+i,T - B'N,j,k+i,r '^^^^ 

(^/yfj.fc+i.rWjv"^^') {'r,yk\y'k) = j VnAvj - yk+i)u''^'^^\T,yk,yk+i;y'k,yk+i)dyk+i, 
i.tUn^^^) {T,yk;y'k) / VW,t(2/j- - yk+l)u']^^^\T,yk,yk+l■,y'k,yk+l)dyk+l■ 



B 



NJ,k+l 



and U^''^ (r) is the solution operator to the free equation, that is 



Then there is a C independent of j, k and N such that 

(•T 







-DAf,i,fc+l,TWjv 



L2 



dT < 



where 



Proof. We prove our main auxiliary theorem in Section 6. This theorem establishes the Klainerman- 
Machedon space-time bound for /3 e (O, |] . □ 

Via the above theorem, we infer that every limit point {uC^)} of juj^^} satisfies the space time bound 

tan uiTn 

r 



i?WB,-fe+in('=+i)(T,.;-) 



L2(R3'=xR3'=) 



dT < C'', 



for some C > and all 1 < j ^ fc. 
Step III. Regarding the solution to the infinite hierarchy 4.2, we have the following uniqueness theorem. 

Theorem 6. Let {u^''^} be a solution of the infinite hierarchy 4-^ in [s,T] c [O, *'^"^^» ] subject to zero 
initial data that is 

u^'\s,yk;y'k) = 0,yk, 

and the space time bound 

i-T 



f iiWB,.fe+it.('=+i)(r,.;-) 

J s 



L2(K3fexR3'=) 

for some C > and all 1 ^ j ^ k. Then Vfc, r e [s, T], we have 



dT < 



(4.5) 



i?('=)««(r,-;-) 



0. 



Proof. See Section 5. 



□ 



Since we have shown the space-time bound 4.5 in Step III, we apply the above uniqueness theorem and 
deduce that 

k _ 

u^'\r,yk;y'k) = [] ^(^. %)^(^. 2/^)' (4-6) 



where ^{t, y) solves the 3D NLS 







2 _ 


bo 








lir) 



in R^+i 



(4.7) 



U0,y) = (j)Q. 

Hence the compact sequence l^jv"*! only one limit point. So 

k 

" V n ~^i'^^yo)4>{T,y'j) as iV ^- 00 

in the weak* topology. Since 77/*^^ is an orthogonal projection, the convergence in the weak* topology is 
then equivalent to the convergence in the trace norm. 
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Example 1. At the suggestion of Professor Walter Strauss, we give a brief explanation on why a factor- 
ized state like formula 4-6 is a solution to the Gross-Pitaevskii hierarchy. Consider k = 1, then plugging 
^(r, ?/i)^(r,yJ) into the infinite hierarchy yields 



idr 



-A,, 



1 , 

-Ay, 

2 ^1 



= bo 



4'] {r,yi)HT,y[) <i>{T,yi) 
bo 



2 — 



{r,y[) 



9{r) g{r) 
boBl2 (^('r.2/i)^(T".2/2)^(r,yi)^('r.2/2)) - boBf ^ (^^(r,yi)^(r,y2)^(T,?/i)^(T,?/^)) 



9{r) 



which is 



9{r) 



boBi^2U^^\ 



Step IV. In Step III, we have concluded the convergence 



Urn Tr 

JV->oo 



,(fe) 



Notice that u}^' = '7]v 



u^^\r,yk;y'k)-l[Hr,yj)kr,yj) = 0, Vr e 
rr^Tw^ and the lens transform of u'^'^^ is 



tanwTn 



^^''\t,^k;^'k) = llcp{t,xj)cl>{t,x'j), 



where ^(f, x) solves 



idtj 



1 



A, 



v ^ 

(/)(0,x) = (j)o{x) 
i.e. equation 1.4. Thence we conclude that 

k 



' + bo \(t)f(t> in M^+i 



Urn Tr 



i=i 



= 0, Vte [0,To], 



as a result of the fact that the lens transform preserves the trace norm (Lemma 4). Since equation 
1.2 is time translation invariant, we have established Theorem 4 and consequently Theorem 1. The 
purpose of the rest of this paper is to prove Theorems 6 and 5. 

Remark 1. Since cp = M^^<j), the global well-posedness of equation 4-7 is implied by the global well- 
posedness in the scattering space ^ of equation 1.4 which comes from the Strichartz estimates. Through 
the lens transform, we always have a L"^ solution to equation -^.7. Lemma 6 then shows Vcp G L"^. 

5. The Uniqueness of Hierarchy 4.2 (Proof of Theorem 6) 

In this section, we produce Theorem 6 with a collapsing estimate and the Klainerman-Machedon board 
game. We list these two tools in below. For convenience, we set the coupling constant bo in the infinite 
hierarchy 4.2 to be 1. 

Lemma 7. [30] Assume u^''~^^^ verifies 



1a 1a 
ic'r + TrAy, , , — -Av' 
' 2 yk+i 2 yk+i 



.C^+i) = 0, 



i2(RxM3'=XR3'=) 



then there is a C > 0, independent of j, k, and s.t. 

(i?,-fc+itx('=+i))(r,yfe;y;,) 

Proof. This is Theorem 1.3 of [30]. For some other collapsing estimates, see [10, 12, 22, 28]. 



□ 
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Lemma 8. Assuming zero initial data i.e. 

7.W(s,yfe;y'fe) = 0,Vfc, 

then one can express u'^^^ti, •; •) in the Gross- Pitaevskii hierarchy 4-2 as a sum of at most 4" terms of the 
form 



or in other words, 



,(1) 



(5.1) 



Here Zn+i = (''"2, T3, t„+i), D C [s,ti]", /x^ are a sei of maps from {2,...,n+l} to {l,...,n} satisfying 
Mm (2) = 1 a?"^ Mm(j) < 3 for all j, and 

^(T„+i,M™)(«("+^^(r„+i)) = |^n5(r,+i)j ;7(i)(ri-T2)Bi,2C/('nT2-T3)B^^(3),2... 

C/("H^n - T„+i)i?^„(„+i),„+i(w("+i)(r„+i, •; •)). 

Proof. The RHS of formula 5.1 is in fact an application of Duhamel's principle involving only the inhomo- 
geneous terms since we have zero initial data. The parameter n is the coupling level we take. This lemma 
follows from the proof of Theorem 3.4 in [30]. One needs only notice that factors depending solely on r, e.g. 

1 



commutes with U'^''^ and -Bj,fc+i Vj, k. 



□ 



With Lemmas 7 and 8, we prove Theorem 6. Let = {(ts, ...,t„+i) | {t2,T3, ...,r„+i) G D} where D 
is as in Lemma 8. Given that we have already checked that 

i?«««(so,-) =0, 

applying Lemma 8 to [so,ti] C [s,T] C [0, t:SSJ£loj^ have 



||i?('^nW(Tl,-) 



E 



i,ri]" 



Z,2(IR3xR3) 

c^Z„+i 



L2(R3x 



V( inf g(T)) / i?W[/«(ri-T2)i?i,2C/(''(T2-r3)B 

i?(l)Bi,2C/('n^2 - T3)B^^(3),2... 



A«™(3),2• 



L2(R3xR3) 



i2(R3xI 



= 

m •'[so.-rij" 

< ^C"(ri-So)^/ (/ i?«Bl,2f/(')(T2-T3)B^^(3),2- 

^ Y,C-CiT^-So)^ f i?(2)B^^(3),C/(3)(^3_^,)... 
TO J[so,ri]"-i 

(Iterate n — 2 times) 



rfr2 1 dT^.-.dr. 



3...aT n+l 



L2(R3xR3) 



L2(R3xR3) 

dT3...dTn+i (Lemma 7) 



™ J Sn 



Z,2(R3xR3) 



dTn 



+ 1 



< C{C{t^-so))— . 
Let (ti — So) be sufficiently small, and n — )• 00, we infer that 



i?«u(i)(ri,-) 



L2(R3xR3) 



in [so,ri] . 
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Such a choice of (ti — sq) works for all of [so,T]. Accordingly, we have ||-R^'°'''"^'^''('''> OH j,2(r3xr3) ~ ^' 
yk,T £ [so,T] by iterating the above argument. Thence we have attained Theorem 6. 

6. The Space-Time Bound of the BBGKY Hierarchy for (3 e (0,2/7] (Proof of Theorem 5) 

We establish Theorem 5 in this section. This section also serves as a simplification and an extension of 
Chen-Pavlovic [7]. Without loss of generality, we may assume k = 1 that is 



J 

Jo 



'■N 



L2 



dT € C. 



(6.1) 



We are going to prove estimate 6.1 for a sufficiently small time T determined by the controlling constant 
in condition 4.4 and independent of A'', then the bootstrapping argument in Section 5 (Proof of Theorem 
6) and condition 4.4 provide the bound for every finite time T e [0, iss^Ia.'^^ Since we work with norms 
here, we transform condition 4.4 into the energy bound: 

2 



n(i-A,j^(i-A,,)^ 



,(fc) 



To obtain the above estimate, one notices 



J {1 - Ay)i {1 - Ay,)i J ^{y,r)^{y',r)dr 
j j {1-Ay)ict>{y,r){l-Ay,)i<f>{y',r)dr 



dykdy', ^ I Tr I [] (1 - A 



dydy' 



Vj) 



dydy' 



j (^j {l-Ayy<k{y,r){l-Ay)"^4>{y,r)dr^ (^j {I - Ay,y 4>{y' ,r) {I - Ay,Y 4>{y' ,r)dr^ dydy' 



i(y,r) (1 - Ay)(j){y,r)dydr 



(k) 

the energy bound then follows from the definition of u}f . The analysis of Theorem 5 also involves i?jv,j,fe+i 
which approximates -Bj^fe+i for every r, so we generalize Lemma 7 to the following collapsing estimate. 

Theorem 7. Suppose ^(t, 1/1,1/2,2/2) solves the Schrddinger equation 

m R^+i 

w(o, 2/1, 2/2, 2/2) = f{yi,y2,y'2), 

and g{T) > cq > 0, then there is a C independent of N and u such that 



IUt 



1a 1a 1a 

-Ay^U+-Ay^U--Ay:^U 



/ \^y \ ( / {9^{r)VN{g{T) (2/ - 2/2))) 5{y2 - 2/^)w(t,2/,2/2;2/2)'^2/2c!2/2 



dydr 



< Chl\\\yyA\Wy,\ 



where bo = JVNdy. 

Proof. Theorem 7 follows from a slightly modified version of the proof of Theorem 2 of [12] . We include it 
in Appendix II for completeness. □ 

We now present the proof of estimate 6.1. In order to more conveniently apply the Klainerman-Machedon 
board game, let us start by rewriting hierarchy 4.3 as 



where 



^ fy{k) (k) , s 

Jo 9\Tk+l) 

N Jo 9{rk+i) 

^irtr"i^^(^'yfc;yfc) = {^N,r{yi-yj)-VN,r{yi-yj))u'k\T-,yk;y'k) 



(6.2) 



(fe+1) _ 



= ^ BNJ,k+l,TUN~^^^ = X] {^N,j,k+1,T 



-B 



N,j,k+1 



(fe+i) 

'at 
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We omit the i in front of the potential term and the interaction term so that we do not need to keep track 
of its exact power. 

Iterate Duhamel's principle (equation 6.2), we have 



= U'- >[T2)Uyo+ U'- '{T2-T3) '— r drsH — U^>{T2-T2,) —, — ^ dT3 



N 

,(3) 



r(3) „(3)/ X 



= '{T2)V>Nfi + — r— / >{t2 - T3) -dry, + \T2 - 

N Jo 9 (T3) Jo g (T4) 



9(^3) 

f,(2) (2) . X 
5(^3) 



+- 



N-2N-3 n 



N N 



Jo 9{t3) Jo 



5(-r4) 

) — ^ dTidT3 



= FreePart'^''^ + Potential Part'^''^ + InteractionPart^'^^ 
where 

FreePart'-''^ 



k / j 



U^'Hr2)u^Z + E[U 

j=3 \l=3 

fT2 fTj-l 



N 



"wo) dT3...dTj, 



X 

PotentialPart^''^ 



= rU^^\T2-T3) 

Jo 



k / j 



9{t2.) 



c^rs+E n 



j=3 V=3 



N + 1-1 
TV 



0) 



V, 



9{Tj+l) 



dTj-^l dT3...dT 



InteractionPart^''^ 



/k+l 

n 

\l=3 



N + l-l 
N 



X (^U^+^'(Tfe+l)) dT3...drfe+i 

Prom here on out, the k's in the formulas are the coupling level we take to prove estimate 6.1, it is distinct 
from the k in the statement of Theorem 5. 
We are going to argue 

dT2 ^ C (6.3) 



/ R^^^BN,i,2,T2FreePart^''\T2) 
Jo 

fT 

/ R^^^BN,i,2,T2PotentialPart^''\T2) 
Jo 

i R''^^ BN,i,2,T2lnteractionPart^''^ (T2) 
Jo 



L2 



1,2 



dT2 ^ C 



dT2 < C 



(6.4) 
(6.5) 
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for some C and a sufficiently small T determined by the controlling constant in condition 4.4 and independent 
of N. We observe that B^N]rj 2j terms inside so that each summand of u^^(r2) contains factorially many 
terms A;!). So we use the Klainerman-Machedon board game to reduce the number of terms. Define 



\l=2 



where Zj^i means (ra, ...,Tj+i) , then the Klainerman-Machedon board game implies the lemma. 
Lemma 9. [30] One can express 

/ ... / JN{Tj+^){f)dTj+^ 

Jo Jo 

as a sum of at most 4-'"^ terms of the form 



/ Jjv(Tj-+i,/i„)(/)(iTj+i, 

Jd 



or in other words, 



[ ... [ Jiv(Tj+i)(/)rfrj+i = 51 / JN{Tj+i,Hm)U)dTj+-^. 

Jo Jo „ Jd 

Here D C [0,T2]-'~^ , jJ-^ are a set of maps from {3, j + 1} to {2, satisfying i2^(3) = 2 and <l 
for all I, and 

,/n^(4),4,r4"- 



• J - ' J+i/-"iv,Ai„(i+l),i+l,r,+l(/)■ 
Remark 2. There is no difference between Lemma 9 and the one we used for the uniqueness of hierarchy 



4.2 (Lemma 8). We have restated it to remind the reader of its exact form since we start from here. 



With the above lemma and the collapsing estimate (Theorem 7), we have the following relation, which 
is essentially part of the proof of Theorem 6, to help establishing estimates 6.3, 6.4, and 6.5. 



/ R^'^^Bn,1.2.T2 / JN{Tj_f_i,II^){f)dTj_^_-^ 

Jq Jd 



dro 



(6.6) 



L2 



I j^\tl9iTl+l)j i?(l)^jv,l,2,r.t/(2)(T2-r3)BjV,2,3,T3-rfr3-dT, + i 

/ X -7+1 

/ i?(')Bjv,l,2,r.C/^''(T2 - T3)Sjv,2,3,T3- 

J[0,T]i 



dT2 



inf (7(r) 



dT2dT3...dTj+i 



/ (/ i?WSw,l,2,r.[/(2)(r2-T3)^w,2,3,r3- 

JfO.Tl'-i \J 



/[O.T] 

(Cauchy-Schwarz) 



L2 



dT2 ] dT3...dTj+i 



< &-^(cT^\j RP'^BN,2,3,r^U^^\Ts-Ti)... dTs...dTj+i {Theorem!) 



'[0,T] 

(Iterate j — 2 times) 



,■1 

{CT^y-^ 

Jo 



3 + 1 f 



1,2 



dn. 



j+i- 



We show estimate 6.3 in Section 6.0.1. Assuming Proposition 3, whose proof is postponed to Section 
6.0.4, we derive estimate 6.4 in Section 6.0.2. Finally, by taking the coupling level k to be InA^, wc check 
estimate 6.5 in Section 6.0.3. We remark that the proof of estimates 6.3 and 6.4 is independent of the choice 
of the coupling level k. Taking the coupling level k to be In AT in the estimate of the interaction part is 
exactly the place where wc follow the original idea of Chen and Pavlovic [7] . Moreover, the proof of estimate 
6.4 (Section 6.0.2) is the only place which relies on /? G (O, |] in this paper. 
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(2) (2) 

6.0.1. Estimate of the Free Part of u}^' . Applying Lemma 9 and relation 6.6 to the free part of m}^ , we 
obtain 

/ R^^^BN,i,2,T^FreePart^''\T2) dra 

^ j=3 m Jo Jd i2 



< CT^ 



-n, Q 



1,2 



1,2 



j=3 m 

k 



j=3 



1,2 



+ C^(CT5)J-iC^ (Condition 4.4) 



j=3 



< C < 00 for T small enough. 
Whence, we have shown estimate 6.3. 

6.0.2. Estimate of the Potential Part of . The same procedure in Section 6.0.1 deduces 



[ Bjv,i ,2,r2 PotentialPart'^''^ (t2 ) 

Jo 



1,2 



dro 



[ / 'i?(^)SAr,i,2,x.C/(2)(T2-T3)- 

Jo Jo 



5(^3) 



drr- 



dro 



L2 



j=3 m 



dTj+i)dT, 



dT-2 



1,2 



jlY.c{CT-y-' f 

A-r, ^ Jo 



j=3 m 



Jo 



n(2)T>(2) „(2), ^ 



1,2 



1,2 



L2 



1,2 



Assume for the moment that we have the estimate 

i-T 



Jo 



1,2 



dT ^ CqC^T 



where C and Co are independent of T, k and A?', then 

i-T 



I R^^^BN^i,2.T2PotentialPart'^^\T2) 
Jo 



L2 



dT2 ^ C < 00. 



for a sufficiently small T and a C independent of k and N. As a result, we complete the proof of estimate 
6.4 with the following proposition. 

Proposition 3. Assume (3 G (Oj f ] cmd V is a nonnegative L^(M^) Ci H^{&.^) function. Then, given 
T e [0, i2iL^:ZL.j ^ ^/lere are C and Co independent ofT, k and N such that 



1,2 



dT ^ CqC^T. 



I R^'^'inPir) 
Jo 

Proof. The proof is elementary and we relegate it to Section 6.0.4. 

Remark 3. This proposition is exactly the reason we restrict (3 G (O, |] in this paper. 



□ 
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(2) 

6.0.3. Estimate of the Interaction Part of u}f' . We proceed like Sections 6.0.1 and 6.0.2. 

/ \\R^^^BN,i,2,T2liT'teractionPart^''\T2) dT2 
Jo ' ' ' 

m Jo Jd 
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dT2 

dr k+i- 



Then the next step is to investigate 



I R'^'^^Bn^i^ (^k+l),k+l,Tk+iUN~^ (Tfe+l) dTk+1. 

Jo ^ 
Without loss of generality, set n^{k + 1) = 1 and look at -Bjv,i,fc+i,rfe+i) ^® have 



J \R^''^B],^i,k+i,r,^yN^'\rk+i)fdykdy', 

= j R'-''^ j VN,Tk+i{yi-yk+i)u''^~^^\Tk+i,yk,yk+i;y'k,yk+i)dyk+i 



^ c 



dykdy'k 

, Yk, j/fc+i; y'fe, yk+i)dyk+i 



dykdy'k 



+c J j VN,rk+i{yi-yk+i)R^''^u''^~^^\Tk+i,yk,yk+i;y'k,yk+i)dyk+i 

= C{I + II). 



dykdy'k 



Noticing that 
/ 



Vr 



Jf2 



^ C ||V"||^2 uniformly for r € [O, iss^] ^ can then estimate 



yN,T^^,\ ivi-Vk+i) Yl \^y^ n l^y'j u^N^^\Tk+i,yk,yk+i;y'k,yk+i)dyk+i 



U=2 



dykdy'k 



^ / dykdy'k 



(VAT.Tfc+i) iVl - Vk+l) 



dyk+1 



k 



n l^f^l n V^y'i {Tk+i,yk,yk+i;y'k,yk+i) 



U=2 



k 



dyk+1 (Cauchy-Schwarz) 



< CN^^ \\V'\\\-, 

< CN'>P\\V'tj^. j dykdy'k 



n l^wl ) I n \^y'o 1 ^N'^^\'rk+i,yk,yk+i;y'k,yk+i 



j=2 



dyk+1 



dykdy'k 



I 



(1 - ^m+.f (i - ^vu)^ n l^^.l n ^y'\\ u^N'^^\Tk+uyk,yk+i\y'k,y'k+i) 



U=2 



dyk+idy'k 



+1 



(Trace Theorem) 
< CiV^/^ II V' II ^2 C'^+i (Condition 4.4) 

and 



= j j '^N,Tk+i{yi-yk+i)R^^^u-N''^\Tk+i,yk,yk+i;y'k,yk+i)dyk+i 

< ||y||2^ (sj^jj^g method as /). 



dykdy'k 
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Accordingly, 



Remark 4. The estimates of I and II may not be optimal. But they are good enough for proving estimate 
6.5 for arbitrary /9 > 0. 



Thence 



dro 



,ik+l] 



(Tk+l) 



dTk+i- 



I \\R^'^^BN,i,2.TjnteractionPart'^''\T2) 
Jo II 

< ^^-^C{CT^f-^T^ (ciV^ \\V\\h, C'^+i) 
C\\V\\H2{T'^fN'^C^. 
Take the couphng level fc = InA^, we have 

/ R^^'^BN,i,2.TjnteractionPart'<''\T2) dT2 \\Vfu^ [T^""^ N^-^ N". 
Jo ' ' 

Selecting T such that 



ensures that 
and thence 



R^^^BN,i,2.r2 InteractionPart^^'' (t2 ) 



L2 



dro < C 



where C is independent of A''. 

We remind the reader that, at this point, wc have obtained estimates 6.3, 6.4, and 6.5 for a sufficiently 
small T determined by the controlling constant in condition 4.4 and independent of A''. Thus one can 
repeat the argument to acquire the estimates for any finite time T G [O, ] through bootstrapping and 

condition 4.4. Whence, we have earned estimate 6.1 and established Theorem 5. The rest of this section is 
the proof of Proposition 3. 

6.0.4. Proof of Proposition 3 . We will utilize the lemma. 
Lemma 10. [18] 



x^-X2)\f{x^,X2)Vdx^dx2^C\\V\\j^, / (1-A,J= (l-AxJV(a^i,a^2) 



In particular, 



j \V{XI - X2)\' \f{xu X2f dXidX2 < C \\V\\l, 



(l-A,JMl-Ax,)V(a;i,a;2) 



dxidx2- 



dx\dx2. 



Proof. This is Lemma A. 3 in [18]. 

Without loss of generality, we show Proposition 3 for A; = 2 which corresponds to 



(I' 



and 



V2 1 \ ^y{ I 



J/2 



V-2 



VN,T{yi-y2) {2), ,x 
MV(^'y2;y2) 



VNAyi-y'2) (2)/ 1^ 

wV('r,y2;y2) 



dT ^ CnC^T 



□ 



(6.7) 



i2 



dT < CqC^T. 



By similarity we only prove estimate 6.7. For a general fc, there arc 2fc2 terms in V^^, whence Rn^VnWn 
has about 8fc^ terms by Leibniz's rule. Since Sk"^ can be absorbed into C*^, the method here applies. 
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First of all, /3 S (O, |] implies the following properties of Vn,t/N '■ 

Vr , 

= iV-(^-H(4-f>o) 



21 



VN,r/N 

/N 

(yN,r)" IN 



= Af3/3-l 



Vr 



< iV-7 



{V, 



LP 



L2 



Vr 



LP 



Vr 



, decays up to p = 6, 



LP 



L^ 



where eq is (| — /3) ^0. On the one hand, 



Vr 



< C \\V\\h2 uniformly for r € [O, is!!^] . On the other 



hand we assume F is a nonnegative L (R'^) n H (W^) function. Thus we know Vr € L°° and 
We compute 



(|VyJ|Vy,| 



^1 1 



= (|V,J|V,J7V-iVW,.(2/i-y2)) (|V,j| 



(iV-'t^iv,x(yi-2/2)w^^(T,y2;y^)) 
V 



2/2 



+( 




N- 


'VnAvi 


-y2)) 


(|V,.| 






+( 




N- 


^VnAvi 


-2/2)) 


(|V,J 


|v.il 





+iV"'V;V,r(2/l - 2/2) (|V„J |Vj,,| |V,,;| 



V2 



(2),' '\ 



Mi^^(T,y2;y2)) 
Wj^^(-r,y2;y2)) 

(2)/- /\ 

wV('^'y2;y2)- 



But 



(|V,J|V,,|Ar-ly;V,.(yi-?/2)) (|V,; 
= sup j (VAr,r)" (2/1 - 2/2) 

(vv,x)"(y) 



2 






|V.;| 



2/2 



Mj^^(r,y2;y2) dy2rfy2 



(-r,y2;y2)) 



2 

L2 



^ sup 

r 

(Lemma 10) 



dy 



(l-A^JMl-A,,)^ V 



wy(T,y2;y2) t^y2 ) rfya 
< C7V-^^° imiL sup 1 |(l-A,JMl-A„J^(l-Aj,;)^(l-A^,)'ti^)(r,y2;y^)|'dy2rfy^ 



< Ar-7eo \\v"\\l2 (Condition 4.4) 
and 



sup 



iV-'t>^,.(2/i-2/2)(|V,J|V,j|V,;| 



y2 



(2 / /\ 

■wV('r:y2;y2) 



= sup j |A^-'VAr,r(yi -2/2)|' |(|V^J |Vy,| IV^^] |v^^|) u^^(r,y2;y^)|' dy2c;y^ 



< sup 

r 



^00 sup 



(1 - A,J^ (1 - A,,)^ (1 - A,;)^ (1 - Ay,y u^^\T,y2;y'2)\ dy^dy'^ 
y |(1 - A,JM1 - A,J^ (1 - A^;)^ (1 - A^,)^ «(^)(T,y2;y^)|' dy2dy^ 



U2 (Sobolev and Condition 4.4). 



Since the same method applies to the middle terms, we have obtained estimate 6.7 and hence Proposition 
3. 

7. Conclusion 

In this paper, we have rigorously derived the 3D cubic nonlinear Schrodinger equation with a quadratic 
trap from the TV-body linear Schrodinger equation. The main novelty is that we allow a quadratic trap in 
our analysis and the main technical improvements are the simplified proof of the Klainerman-Machedon 
space-time bound as compared to non-trap case in Chen and Pavlovic [7] , and the extension of the range 
of /3 from (0, 1/4) in Chen and Pavlovic [7] to (0, 2/7] . Compared to the 2D work [12] which is also by the 
author, the 3D problem in this paper is of critical regularity. To explain what we mean by critical: in 2D 
one easily obtains the |V| ^-space-time bound needed for the uniqueness theorem by a trace theorem: in 3D 
the only way to obtain the space-time bound 4.5 is through smoothing estimates since one does not have 
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enough regularity to apply a trace theorem. Thence the key arguments in 3D are more involved and totally 
different from the 2D case which is a subcritical problem. Moreover, we have established the trace norm 
convergence in the main theorem which is a stronger result than the Hilbert-Schmidt norm convergence. 
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9. Appendix I: Proof of Corollary 1 

For the purpose of this Appdendix I, we may assume w = 0. Or in other words, we skip Steps I and 
IV of the proof of Theorem 4 here. When the desired limit is an orthogonal projection, one does not need 
this proof. This is a functional analysis argument and all operators mentioned in this Appendix I acts on 
jji ^^k\ ^ -yyg prove CoroUary 1 by verifying the hypothesis of the following lemma. 



Lemma 11. [35] Assume the operator sequence {An} satisfies that, as bounded operators, An A, An 
A* and ^ \A\ in the weak sense. If 



then 



lim Tr|A„| =Tr|A| 

JV— >oo 



lim Tr|^„ - A] = 0. 



Proof. This is Theorem 2.20 in [35]. It implies the Griimm's convergence theorem (Theorem 2.19 of [35]) 
used in [21]. □ 

We first observe that condition 1.6 implies the a-priori estimate 

k 



sup Tr m(l-^-.) hiv'^^C-^ 



te[o,T] 



Thus we have the compactness argument and the uniqueness argument to conclude that, as trace class 
operator kernels, 

7^^(i,Xfe;x;) - 7W(i,Xfe;x;) (weak*). 

Remark 5. Because we assume w = 0, both of the Erdds-Schlein-Yau uniqueness theorem [18] and the 
Klainerman- Machedon uniqueness theorem [30] apply here. For the general case, one has to use the main 
argument in this paper. 

The above weak* convergence as trace class operator kernels infers that as Hilbert-Schmidt kernels, 

7i^^(t,xfe;x',) -7W(t,Xfe;x',) (weak*), 

because there are less test functions. Since L'^ {dx^dx'/.) is reflexive, the weak* convergence is no different 
from the weak convergence. Thus we know that as Hilbert-Schmidt kernels and hence as bounded operator 
kernels, 

7^'^(t,Xfe;x'J -7<')(txfc;x'fc) (weak). 
At this point, we have verified that, as bounded operators. An — ^ A and A* ^ A* in the weak sense since 
7^^ and 7^*^^ are self adjoint. We now check — ^ |^| . 

To check \An\ — ^ j^j, one first notice that 7^^^ and 7^^°^ has only real eigenvalues since they are self 
adjoint. Moreover, 7j^^(f, Xfc;x^) has no negative eigenvalues by definition 1.3, in fact, 

2 



/(/ 



^{x, z)4>{y, z)dz^ f{y)dy^ f{x)dx = j dz 



(a;, z)f{x)dx 



>0. 
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Notice that /(xfe)/(x^) in the estimate above is also a Hilbert-Schmidt kernel, the definition of weak 
convergence in L^(dxferfx^) then gives 

J 7^'^ni,Xfe;x;)/(xfc)/(x'Jdx,dx; = ^irn^J {t,^k;^'k)Wk)fK)d^kd^'k > 0, 

as (i, Xfe; x'^,) — ^ 7^'^^ {t, x^; x^) weakly in L^(dxfedx^). So we have checked |A„| — ^ \A\ because |A„| = An 
and \A\ = A. 

To prove Corollary 1, by Lemma 11, it remains to argue 

7^^(t,Xfe;x;) 



lim Tr 

N->oo 

Notice that we have the conservation of trace 



Tr 



j -fP{t,Xk;Xk)dxk = y 7^' (0,Xfe;Xfe)c;xfc (By definition 1.3) 
j 7('=)(t,xfe;xfe)(ixfc = j 7W(0,Xfe;Xfe)dxfc [8]. 
and we have shown that 7^^(t,Xfe;x'j,) and 7('^)(t, x^; x'^,) have no negative eigenvalues, hence 
7jv^(i,Xfe;x^) = y"7j^^(t,Xfc;Xfc)c;xfc = y 7^^(0,Xfe;Xfe)dxfc, 

7('=)(i,Xfe;x;) = y7W(t,Xfc;Xfe)dxfc =y7('=)(0,xfe;xfe)rfxfe. 
On the one hand, 

y 7('H0, Xfe; Xfc)dxfe = Tr I7W (0, Xfc; x',) 
on the other hand. 



Tr 
Tr 



lim Tr 

JV->oo 



7^ (0,Xfe;x'fe) 



lim Tr 

JV->-oo 



Tr 



v(fc) 



(*,Xfe;Xfe) = y 7('')(f,Xfc;Xfc)dxfc = j 7^'') (0, x^; Xfe)dxfc. 



That is 



Tr 



7<')(i,x,;x;,) 



lim Tr 

Af-s-oo 



7^^(i,Xfe;x^) 



Whence we conclude the proof of Corollary 1 by Lemma 11. 

10. Appendix IT Proof of Theorem 7 

In this Appendix II, instead of proving Theorem 7, we establish the following general estimate. On the 
one hand, the special case (Theorem 7) does not yield a shorter proof. Unlike estimate A. 8 in [7], Theorem 
7 does not follow directly from [29, 30]. The time dependence of the potential makes the proof closer to the 
Beals-Bezard simplification [2] of the Klainerman-Machedon null form paper [29] . On the other hand, the 
anisotropic version of the main theorem (Theorem 3) requires the following theorem. 

Theorem 8. Define the diagonal mMrices 

Ui{t) 






where ai G Ljg^{M.) functions and 



a2{t) 
asit). 



and Bt 




ai,f5i ^ Co > a.e. 



Suppose M(i, xi, X2, Xg) solves the Schrddinger equation 

iut + divxi {LtVxJu + diva;2 (itV^^s) u - div^^^ (^tV^^) u = in 



(10.1) 



u{0,Xi,X2,x'2) = f{xi,X2,X2), 



then there is aC independent of N and u such that 

IV 



/ |v.|(/ 

'0 



< Chf 



where 



Vlf{x — X2)5{X2 — x'2)u(t, X, X2'i X2)dX2dX2 
I I I /II2 ) 

V^{x)=det {Bt)VN{Btx), 



dxdt 
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which has the property that 



sup 



< / \V^{x)\dx = bo. 



Through out this paper, we have used x to denote a 3D variable. However, the anisotropic feature of 
Theorem 8 requires us to use the individual components of a 3D vector in the proof. So we use bold symbols 
to represent a 3D vector in this Appendix II. e.g. x = {xi,X2,X3)gM.^. 

We will make use of the lemma. 

Lemma 12. [30] Let ^ e and P be a 2d plane or sphere in with the usual induced surface measure 
dS. 

(1) Suppose < a, 6 < 2, a + 6 > 2, then 



dSiv) 



c 



.\a+b-2 ■ 



(2) Suppose e = then 





1* |2-£ 1 |2-£ 


l-v 





C 



Both constants in the above estimates are independent of P. 
Proof See pages 174 - 175 of [30]. 



By duality, to gain Theorem 8, it suffices to prove that 

j h{t,x) |Vx| V^(x- X2)^(x2-X2)u(i,x,X2;x2)dx2dx2^ dxdt 

< C6o||/l|l2||VxiVx,Vx^/||2. 

For this purpose, let 

/j*aiis)ds ^ 

At= \ /o a2(s)ds 

V Joa3{s)ds) 

Consequently, the solution of equation 10.1 is recast to 

n(t,xi,X2,x^) = I e^(«"^'«i+«"^'«=±(^^)"^*«^)e-^«^e-=«=e-^«^/(^i,^2,^2)rf^irf^2rf€2, 
which implies 

j V^{X - X2)(5(x2-X2)w(i,X,X2;X2)(iX2(iX2 

= / ^(^2)e^^^^^'^^+^^^*^^^(^^)"^*^^^e-(«^+«^+«^)/(^i,^2,^2)rf^irf^2rf^2- 
Accordingly, the spatial Fourier transform of 

|Vx| (^J V;(r(x- X2)(5(x2-X2)u(t,X,X2;X2)dX2dX2 



IS 



ICll / ^^(^2)e^^^^^-^^-^^)"^'^^^-^=-^^^+^"^'^=^(^^^"^'^^^/(^l -^2-€2>^2,l2)'^^2rf€2. 



□ 
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which aUows us to compute 



J h{t,x)\Vx\ (^J V^(x - X2)^(x2-X2)u(i,X,X2;X2)rfX2dx'2^ dxdt 



(spatial Fourier transform on h) 



^ /(/i) ||VxiVx,Vx^/||^, (Cauchy-Schwarz) 



where 



-di-^di2dS,'2- 



1^1 ~ ^2 ~ ^2! 1^2! 1^2! 

So our purpose in the remainder of this Appendix II is to show that 

Noticing that, away from the factor V)^(^2)' the integral I{h) is symmetric in ~ ^2 ~ ^2! ^^^^ |^2l > 

suffices that wc deal with the region: |^]^ ^ ^2 ^ ^2! > 1^21 only since our proof treats V^(^2) ^ ^ harmless 
factor. We separate this region into two parts, Cases I and II. 

Away from the region |^i — ^2 "^2] > l^2l) there are other restrictions on the integration regions in 
Cases I and II. We state the restrictions in the beginning of both Cases I and II. Due to the limited space 
near "/", we omit the actual region. The alert reader should bear this mind. 

10.0.5. Case I: I{h) restricted to the region |^2| < 1^21 integration order d^2 prior to d^^- Write the 
phase function of the dt integral inside I{h) as 



(^1 - ^2 - ^2) A, (^1 - ^2 - Q+ilAti2 ± (€2) ^tl2 



The change of variable 
leads to the expression 

m 



^2,new *^2,old 



^1-^2 



(10.2) 



/ 



/V^A(^2)e' 



S2 2 



?2 H 2 



<id^2d^2 



S2 2 



S2 + 2 



1^2 



N 



^^^^^,(2(iizlir^£(iiziii+2€^A,€,±(«i)^A,«^) 



+2C-A„€,±(«Q ^"«^)/i(t,^i)Fj^'(C2)/i(t',^i)c;i(^i'c;€W^2t^4: 



y d^iy J{h){t,^,)h{t,ii)dt 



where 



S2 2 



t 1 €1-^2 
S2 ^ 2 



1^2 



A 2 ±(€2) (^'-^t')€2)/j(i/^|^ 



)dt'd^2d^'2. 
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Assume for the moment that 



/ 



dt ^ Chi 



2 



with C independent of /i or then we deduce that 

Hence we end Case I by this proposition. 
Proposition 4. 

j \J{f){t,i,)\'dt<.Cbl \\f{;i,)\\l. 

where C is independent of f or 



Remark 6. To avoid confusing notation in the proof of the proposition, we use f{t', ^i) to replace h{t', ^i). 
Proof. Again, by duahty, we just need to prove 



J{f){t,^,)g{t)dt 



^c'||/(.,Ci)ILj||5lL- 



For convenience, let 



Then 



J J{f){t,i,)g{t)dt 



/ 



?2 2 



€2 + ^ 



1^21 



{vl{^2)e-''^'-''^''^''^9{t))dtdt'd^2d^'. 



S2 2 — 



S2 H 2 — 



Illl'rfg2 

i€2r 



dl2 



t ii is. 



^2 H 2~ 



To deal with the dt and rft' integrals, for every fixed let 



then 



du „ai (0^2,1 + a2(t)?2,2 + ^3(0^2,3 ^ „ ^„ 
^=2 >2co>0 



which provides a well-defined inverse t{u). 
Consequently, the integral 



1^2 



dt 



is indeed the Fourier transform of 



dt 
du 



du 



dt 
du 
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This is well-defined since 



du = 



vli.i2)'w)\ 



g-i</.(t(u),€i,$2)g(t(u)) 

dt 



dt 
du 



du 



dt 
du 



< sup 



vh{i2) 1 / \9{t)Y 



dt 
du 



dt 



Hence, we find 

/ j{f){t,^,)mdt 



S2 2 



S2 + 2 



1^2 





G(|^2l'mi^2l%^l) 




S2 2 


2 


S2 + 2 


2 



M^sup 









t €i-«2 

?2 2 


2 


€2 + ^' 



p^dpdcr (spherical coordinate in ^2) 



/ 



1^2 



p^da 



V p 



?2 2 



t J. ^l~^2 

S2 H 2 



pdp 



G{p') 



pdp 



(Holder in p) 
< Cbl\\f{;^,)h.\\9h.{ 

However, we have the estimate 



1^1 



sup 



^2 P 



p^dcr 



nil 



2 

2 sup 



^2 " 



C _ ^1~^2 

S2 2 



(?da 



t ^l~^2 

S2 2 



tf _l_ ^l~^2 

?2 + 2 





^2- 


2 


2 


t ^i-€2 

?2 2 




^1-^2-^2 


Is2i 



(Reverse the change of variable in formula 10.2 



^2 + ^ 



\ 

! 



\ 1 

r c?l2 



rf^2 



V 

^^2 





e €i-€2 
^2 2 


2 


t ^1-^2 

?2 2 


IC1-C2- 


-^^r-^ic2r 



I 



2+2e 



/ |3-2e 



(Second part of Lemma 12) 



In the above calculation, the <t in the first line lies on the unit sphere centered at the origin while the cr 
in the second line is on a unit sphere centered at ^^^^ . We use the same symbol because Lebesgue measure 
is translation invariant. 
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Thus, we conclude that 



/ j{m,^,)9it)dt 



□ 



10.0.6. Case II: I{h) restricted to the region |^2| > l^2l with integration order d^2 P^ior to d^2- Fo^' this 
case, we express the phase function as 

(^1 - ^2 - ^2^^ (^1 - ^2 - ^2) + i2M2 - {i'2fM'2 

= (^1 - ^2)^ At (Ci - ^2) - 2 (^1 - ^2)^ M2 + ilM2 



and let 



J 



{h){t,^,) 



|4i I ^ e-2'(«i -«2)'^^t«2 e2'(€i -^2)^ A'€2 



^l-C2-r2r 1^21' 1^2 



2 1^ |2 1^/ |2 
2 



^i^(€2)^^A'(^2)e-'^(*''^-^=^e-'^(*.«i.€i)/i(t',^i)di'd^^d^2- 
Again, wo want to prove the following property. 

Proposition 5. 

/ \Jm,^l)\'dt^Cbl \\fi;^,)\\% 
where the constant C is independent of f or 
Proof. We calculate 

J{f){t,$r)W)dt 

||^|2g-2j(€i-C2)^At€2e2*(€i-«2)^^t'C2 



2 1^ |2 1^/ |2 
2 



1^1 -^2-|2r 1^21' 1^2 

(^^(42)e-'*(*'«-«^)5W) dtdt'd^^di: 



di2 



|4l ^2 ^2! 1^2 1 



y e2i(€i-«2)"A*'€2 (|l/j^'(^2)e-'^^*''^^'^=V(t',^i)) dt' 

Fix — ^2 ^-nd ^2, and write 

where u) = {uj-i,uj2,ujs) is a unit vector in R^. Without loss of generality, we assume that 

max{|a;i| , \uj2\ , jwal} = |wi| 



which in turn implies 



We then write 



1 

V3 



^ \uJi\ ^ 1. 



$'2 = (a;,0,0) + (0,2/1,^2) 

i{t) = 2 / ai{s)ds. 
Jo 



Again u is invertible with 



du 



^ 2co > 0. 
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So we have 



-m(c^i|«i-€2|,T) ^f^,-^^)g-2*|«i-€2l(O,a;2,a;3)^^t(„)(O,yi,y2)g-*0(f(«),€i,«^)^(-^(-y)) ^ ^ 



= c;(-u;i|^,-^2|a;) 



where 



The last expression still has the property that 

J \G{u)fdu^^J\g{t)fdt. 
Just as in case 1, this procedure furnishes 

/ Am,ii)W)dt 



dt 
du 



d^. 



d^'2 



4l ^2 ^2! 1^2 



j e-2i(€i-€2)^A,€i (l^A(^2)e-*'^(*'«i'«=)5(i)) 



dt 



dxdyidy2 



1^1 ^2 ^2! 1^2! 

dxdyidy2 
|Ci ~ ^2 ~ ^2! 1^2! 

sup 



72 ^^2 



dy\dy2 



I '^1 1 1^1-^211^: 



72 ^^^2 



1^1-^211^2 



1^1 ~ C2I 1^2 



dx 



G{x) dx I d$,2 (Holder in x) 



l^il 



sup 



dyidy2 



ICl ^2 ^2! 1^2 



rf^2- 



2 1^1 -^211^2! 

The first part of Lemma 12 along with the restrictions |Ci — C2 — $2! > IC2I IC2I < IC2I entail 

|2 / 



2 1^1-^21 1^2! 

2 



sup 



dyidy2 



^2 ^2! |^2| 



< c 

which finishes the proposition. 



Ilil 



2 1^1-^21 1^2 
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1^1 -^2-e 



/ |2-£ I ^/ l2-£ 
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2 1^1- ^2 r'' 1^2 r'' 



References 



□ 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Observation of Bose-Einstein 

Condensation in a Dilute Atomic Vapor, Science 269 (1995), 198-201. 
[2] M. Beals and M. Bezard, Nonlinear Field Equations: Not Necessarily Bounded Solutions, Journees equations aux derivees 

particUcs 20 (1992), 1-1.3. 

[3] R. Carles, Nonlinear Schrodinger Equation with Time Dependent Potential, Commun. Math. Sci. 9 (2011), 937-964. 
[4] S. L. Cornish, N. R. Clausscn, J. L. Roberts, E. A. Cornell, and C. E. Wieman, Stable ^^Rb Bose-Einstein Condensates 
with Widely Tumable Interactions, Phys. Rev. Lett. 85 (2000), 1795-1798. 



30 



XUWEN CHEN 



[5] T. Chen and N. Pavlovic, On the Cauchy Problem for Focusing and Defocusing Gross- Pitaevskii Hierarchies, Discrete 

Contin. Dyn. Syst. 27 (2010), 715-739. 
[6] T. Chen and N. Pavlovic, The Quintic NLS as the Mean Field Limit of a Boson Gas with Three-Body Interactions, J. 

Punct. Anal. 260 (2011), 959-997. 
[7] T. Chen and N. Pavlovic, Derivation of the cubic NLS and Gross- Pitaevskii hierarchy from manybody dynamics in d = 2,3 

based on spacetime norms, arXiv:1111.6222. 
[8] T. Chen, N. Pavlovic, and N. Tzirakis, Energy Conservation and Blowup of Solutions for Focusing Gross-Pitaevskii 

Hierarchies, Ann. I. H. Poincare 27 (2010), 1271-1290. 
[9] T. Chen, N. Pavlovic, and N. Tzirakis, Multilinear Morawetz identities for the Gross-Pitaevskii hierarchy. Contemporary 

Mathematics 581 (2012), 39-62. 
[10] X. Chen, Classical Proofs Of Kato Type Smoothing Estimates for The Schrddinger Equation with Quadratic Potential in 

]R"+i with Application, Differential and Integral Equations 24 (2011), 209-230. 
[11] X. Chen, Second Order Corrections to Mean Field Evolution for Weakly Interacting Bosons in the Case of Three-body 

Interactions, Arch. Rational Mech. Anal. 203 (2012), 455-497. DOI: 10.1007/s00205-011-0453-8. 
[12] X. Chen, Collapsing Estimates and the Rigorous Derivation of the 2d Cubic Nonlinear Schrddinger Equation with 

Anisotropic Switchable Quadratic Traps, J. Math. Purcs Appl. 98 (2012), 450-478. DOI: 10.1016/j.matpur.2012.02.003. 
[13] P. Cladc, C. Ryu, A. Ramanathan, K. Hclmcrson, and W. D. Phillips, Observation of a 2D Base Gas: From Thermal to 

Quasicondensate to Superfluid, Phys. Rev. Lett. 102 (2009) 170401. 
[14] K. B. Davis, M. -O. Mcwcs, M. R. Andrews, N. J. van Drutcn, D. S. Durfcc, D. M. Kurn, and W. Kettcrlc, Bose-Einstein 

condensation in a gas of sodium atoms, Phys. Rev. Lett. 75 (1995), 3969-3973. 
[15] A. Elgart, L. Erdos, B. Schlein, and H. T. Yau, Gross-Pitaevskii Equation as the Mean Field Limit of Weakly Coupled 

Bosons, Arch. Rational Mech. Anal. 179 (2006), 265-283. 
[16] L. Erdos and H. T. Yau, Derivation of the Non-linear Schrddinger Equation from a Many-body Coulomb System, Adv. 

Theor. Math. Phys. 5 (2001), 1169-1205. 
[17] L. Erdos, B. Schlein, and H. T. Yau, Derivation of the Gross-Pitaevskii Hierarchy for the Dynamics of Bose-Einstein 

Condensate, Comm. Pure Appl. Math. 59 (2006), 1659-1741. 
[18] L. Erdos, B. Schlein, and H. T. Yau, Derivation of the Cubic non-linear Schrddinger Equation from Quantum Dynamics 

of Many-body Systems, Invent. Math. 167 (2007), 515-614. 
[19] L. Erdos, B. Schlein, and H. T. Yau, Rigorous Derivation of the Gross-Pitaevskii Equation, Phys. Rev. Lett. 98 (2007), 

040404. 

[20] L. Erdos, B. Schlein, and H. T. Yau, Rigorous Derivation of the Gross-Pitaevskii Equation with a Large Interaction 

Potential, J. Amer. Math. Soc. 22 (2009), 1099-1156. 
[21] L. Erdos, B. Schlein, and H. T. Yau, Derivation of the Gross-Pitaevskii Equation for the Dynamics of Bose-Einstein 

Condensate, Annals Math. 172 (2010), 291-370. 
[22] M. G. Grillakis and D. Margctis, A Priori Estimates for Many-Body Harmltoman Evolution of Interacting Boson System, 

J. Hypcrb. Diff. Eqs. 5 (2008), 857-883. 
[23] M. G. Grillakis, M. Machedon, and D. Margctis, Second Order Corrections to Mean Field Evolution for Weakly Interacting 

Bosons. I, Commun. Math. Phys. 294 (2010), 273-301. 
[24] M. G. Grillakis, M. Machedon, and D. Margetis, Second Order Corrections to Mean Field Evolution for Weakly Interacting 

Bosons. II, Adv. Math. 228 (2011) 1788-1815. 
[25] E.P. Gross, Structure of a Quantized Vortex in Boson Systems, Nuovo Cimento 20 (1961), 454^466. 
[26] E.P. Gross, Hydrodynamics of a super fluid condenstate, J. Math. Phys. 4 (1963), 195-207. 

[27] W. Ketterle and N. J. van Druten, Evaporative Cooling of Trapped Atoms, Advances In Atomic, Molecular, and Optical 
Physics 37 (1996), 181-236. 

[28] K. Kirkpatrick, B. Schlein and G. Staffilani, Derivation of the Two Dimensional Nonlinear Schrddinger Equation from 

Many Body Quantum Dynamics, Amer. J. Math. 133 (2011), 91-130. 
[29] S. Klaincrman and M. Machedon Space-time estimates for null forms and the local existence theorem. Comm. Pure Appl. 

Math. 46 (1993), 1221-1268. 

[30] S. Klainerman and M. Machedon, On the Uniqueness of Solutions to the Gross-Pitaevskii Hierarchy, Commun. Math. 
Phys. 279 (2008), 169-185. 

[31] E. H. Lieb, R. Seiringer and J. Yngvason, Bosons in a Trap: A Rigorous Derivation of the Gross-Pitaevskii Energy 

Functional, Phys. Rev. A 61 (2000), 043602. 
[32] E. H. Lieb, R. Seiringer, J. P. Solovej and J. Yngvason, The Mathematics of the Bose Gas and Its Condensation, Basel, 

Switzerland: Birkhaiiscr Vcrlag, 2005. 
[33] O. Penrose and L. Onsagcr, Bose-Einstein Condensation and Liquid Helium, Phys. Rev. 104, 576-584 (1956). 
[34] L.P. Pitaevskii, Vortex Lines in an Imperfect Bose Gas, JETP 13 (1961), 451-454. 

[35] B. Simon, Trace Ideals and Their Applications: Second Edition, Mathematical Surveys Monogr. 120, Amer. Math. Soc, 
Providence, RI, 2005. 

[36] H. Spohn, Kinetic Equations from Hamiltonian Dynamics, Rev. Mod. Phys. 52 (1980), 569-615. 

[37] D. M. Stamper-Kurn, M. R. Andrews, A. P. Chikkatur, S. Inouye, H. -J. Miesner, J. Stenger, and W. Ketterle, Optical 

Confinement of a Bose-Einstein Condensate, Phys. Rev. Lett. 80 (1998), 2027-2030. 
[38] S. Thangavelu, Regularity of Twisted Spherical Means and Special Hermite Expansions, Proc. Indian Acad. Sci. (Math. 

Sci.), Vol. 103 (1993), 303-320. 

[39] K. Yajima and G. Zhang, Local Smoothing Property and Strichartz Inequality for Schrddinger Equations with Potentials 
Superquadratic at Infinity, J. Differ. Equations. 202 (2004), 81-110. 

Dep.'Vrtment of M.'\tiiem.'\tics, Brown University, 151 Thayer, Street, Providence, RI 02912 
E-mail address: chenxuwen@math.umd.edu / chenxuwen@math.brown.edu 



